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1 . 


SUMMARY OF ACTIVITIES 


During the past several months we have concentrated on further analysis of our field 
measurements of carbon dioxide and water vapor fluxes. These data were examined in conjunction 
with reflectances obtained from a helicopter mounted Modular Multiband Radiometer (MMR) at our 
grassland study site (#16) during 1987. These measurements are representative of the canopy scale 
(10-100 m) 2 and provide a good basis for investigating the hypotheses/relationships potentially useful 
in remote sensing applications. 

All of the micrometeorological data collected at site 916 during FIFE-89 have been processed 
and fluxes of C0 2 , water vapor and sensible heat have been calculated. Soil C0 2 fluxes were estimated 
using the relationships developed at our site by Norman et al. (1991). Employing these soil C0 2 flux 
values, in conjunction with our micrometeorological measurements, canopy photosynthesis is being 
estimated. 

A biochemical model of leaf photosynthesis was adapted to the prairie vegetation. The 
modeled leaf photosynthetic rates were scaled up to the canopy level. Further tests are underway 
employing our measurements made during FIFE-89. This model and a multiplicative stomatal 
conductance model are also used to calculate canopy conductance. Further evaluation of these 
calculations, using the FIFE-89 data, is in progress. 

Our future plans include completion of these tests of our canopy photosynthesis and 
conductance models against our FIFE-89 measurements. Our FIFE-87 and FIFE-89 data will also 
be used to examine the contribution of this temperate grassland in northeastern Kansas as a source 
and sink for atmospheric C0 2 under a wide range of meteorological and soil moisture conditions. 


2. DETAILS OF WORK PERFORMED 

2.1 Photosynthesis and Stomatal Con ductance Related to Reflectance on the Canopy Scale 

2.1.1 Introduction 

In recent years many researchers ( e.g ., Monteith, 1977; Tucker et al ., 1981; Asrar et al., 1984; 
Sellers, 1985, 1987) have investigated the relationships between linear combinations of reflected 
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visible and near infrared radiation known as vegetation indices, and biophysical quantities, such as 
the leaf area index, biomass, and the amount of photosynthetically active radiation (PAR) absorbed 
by the canopy. Two vegetation indices commonly used are: 

Simple Ratio Vegetation Index ( SR) = RJRy 0 ) 


and 

Normalized Difference Vegetation Index (NDVI) = (R N - RJI(R M + RJ , (2) 

where R N and R* are upwelling radiances or surface reflectances in the near-infrared and visible-red 
spectral regions, respectively. 

Canopy conductance (g c ) and canopy photosynthesis (P c ) can be described using two simple 
multiplicative relationships (see Jarvis, 1976; Sellers, 1987): 

g c = 8 ; (PAR) • W) AD) xe;] (3) 

P c = />; (PAR) • [AT) AD) AW < 4 > 


where g' (PAR) and P* (PAR) are the unstressed values of gt and P c , regulated only by vegetation type, 
health and density and the amount of incoming PAR. The stress factors f(T), f(D), and f(0,) account 
for the effects of temperature (T), ambient vapor pressure deficit (D) and extractable soil water (0,). 
The "stress factors" f(x) vary between unity (optimal conditions) and zero (unfavorable conditions 
shutting down photosynthesis and transpiration). The stress factors operating on g t and P c are assumed 
to be identical and invariant within the canopy for a given set of environmental conditions. In this 
study, canopy level flux data and the leaf stomatal conductance of the dominant grass species 
measured at the site were used to obtain the functional forms of the expressions for the stress factors 
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(see Appendix A). The derived coefficients for these stress factors are independent of the time- 
varying leaf area index. 

The terms g* and P* are dependent on the vegetation canopy and the incident PAR flux. Prior 
analyses (e.g., Sellers, 1987) have indicated that g’ and P* should increase almost linearly with PAR 
and, therefore, dgJdPAR and dP'jdPAR can be treated as conservative surface properties. Put 
simply, a bare soil surface should show no change in conductance or "photosynthesis" with increasing 
PAR and so its associated partial derivatives of g* and P* would be zero. A thickly vegetated surface, 
on the other hand, would exhibit a strong sensitivity to PAR (non-zero partial derivatives of g* and 

p;>. 

The theoretical study of Sellers (1987) proposed that for ideal conditions SR should be near- 
linearly related to the absorbed PAR (APAR) and to the derivatives of g* and P* with respect to PAR: 

SR « APAR, — — , — — (*) 

dPAR dPAR 

The relationship described in Eq. (5) should hold where (a) the canopy consists of spatially uniform 
cover of green leaves, (b) the reflectance of the background and understory in the chlorophyll 
absorption band is approximately equal to that of the overlying vegetation and (c) the wavebands of 
the sensing instrument are configured so that the perceived scattering coefficients of the leaves 
conform to: 


»» = 1 " 



( 6 ) 


where and w. are the leaf scattering coefficients in the near inf rared and PAR wavelength intervals, 
respectively, G(p) is the relative projected area of leaf elements in the direction of cos' 1 n and p is the 
cosine of the zenith angle of the incident radiation flux (see Sellers, 1987). 
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During 1987, the First ISLSCP Field Experiment (FIFE) was conducted in a 15 km x 15 km 
area near Manhattan, Kansas. In this paper we report on an analysis of the surface flux and surface 
reflectance measurements made at a prairie site in the FIFE study area. The fluxes of water vapor 
and carbon dioxide were measured using the micrometeorological eddy correlation technique (see e.g., 
Kaimal, 1975; Kanemasu et al., 1979; Verma, 1990), and the reflectance measurements were made 
with a helicopter-mounted Barnes Modular Multiband Radiometer (MMR). The effective narrow- 
band visible and near-infrared scattering coefficients for Landsat TM bands 3 and 4, which 
correspond to the bands on the helicopter MMR, were calculated to be 0.12-0.18 and 0.95-0.97 
respectively for individual grass leaves from the data of Waiter-Shea (personal communication) 
collected during FIFE. For calculations involving the penetration of broad-band PAR (0.4 - 0.7 pm) 
into the canopy we have assumed 1 «, - 0.2, following the calculation of Sellers (1987). Equation (6) 
is, therefore, approximately satisfied for a spherical leaf angle distribution (G(p)-0.5) and so we 
might expect that Eq. (5) should hold reasonably well for the prairie site being studied here. The 
purpose of this paper is to test the hypothesis represented in this equation with the field measurements 
collected in this study. The data presented here are representative of canopy scale and, therefore, can 
be considered as a test case for investigating the use of Eq. (5) for remote sensing applications. 

2.1.2 Materials and Methods 

2.1.2. 1 Site 

The study was conducted in 1987 at a tallgrass prairie site (39° 03’ N, 96° 32’ W, 445 m above 
m.s.l., FIFE site #16) near Manhattan, Kansas. The soil is predominantly Dwight silty clay loam 
(Typic Natrustolls). The vegetation is dominated by three C 4 grass species: big bluestem (Andropogon 
gerardii ), indiangrass ( Sorghastrum nutans) and switchgrass ( Panicum virgatum). The prairie was 
burned in the spring of 1987. The experimental area had been lightly grazed for several years, but 
was not grazed in 1986 and 1987. 


'Values of derived from the observations of E. A. Waltcr-Shca (personal communication) vary 
from 0.17 to 0.19. 
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2.1.2.2 Micrometeorological Measurements 

Fluxes of carbon dioxide, water vapor, sensible beat and momentum were measured with 
eddy correlation sensors mounted 2.25 m above the ground. The sensors included: a three- 

dimensional sonic anemometer (Kaijo Denki Co., Tokyo, Japan, Model DAT-310), one-dimensional 
sonic anemometers (Kaijo Denki Co., Tokyo, Japan, Model DAT- 110; Campbell Scientific, Logan, 
Utah, Model CA27), fine wire (0.025 mm) thermocouples (Campbell Scientific, Logan, Utah, Model 
127), a Lyman-alpha hygrometer (AI.R., Inc., Boulder, CO, Model LA-1) with a 5 mm path length 
and a rapid response C0 2 sensor with a 200 mm path length (Bingham el al., 1978; Anderson el al., 
1984). Net radiation, PAR, soil heat flux, air temperature and humidity were also measured. Details 
on the instrumentation and procedure can be found in Kim and Verma (1990a & b). The 
configuration of the flux instrumentation indicates that an area extending up to (100 m) 1 upwind of 
the instrumentation would have the strongest (~70%)influence on the measured flux (see e.g. Schuepp 
el al., 1990). 


2. 1.2.3 Plant and Soil Measurements 

Stomatal conductance was measured on selected leaf blades of tillers of big bluestem, 
indiangrass and switchgrass, employing a steady-state porometer. The surface soil water content (0- 
0.1 m) was monitored gravimetrically almost eveiy day and the subsurface water content (0.1-1.4 m) 
was measured with a neutron probe on a weekly basis. Leaf area index was measured generally every 
two weeks. 

2.1.2.4 Reflectance Measurements 

Measurements of surface reflectance were made with a helicopter-mounted MMR from 230 
m above the ground. The MMR has eight spectral channels, however only two channels were used 
in the analysis: 0.63 - 0.68 pm (MMR3 or visible-red) and 0.75 - 0.88 pm (MMR4 or near infrared). 
A second MMR was positioned on the surface to measure reflected radiance from a Barium Sulfate 
reference panel as a measure of irradiance. These reference panel measurements were made at one 
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minute intervals starting thirty minutes before the helicopter missions and ending 20 to 30 minutes 
after the helicopter left the FIFE study area. The reference panel was calibrated using procedures 
outlined in Jackson et al. (1987), and the MMR was calibrated using procedures described in 
Markham et al. (1988), thus allowing determination of absolute radiance from the MMR voltage 
responses. 

A minimum of 25 and a maximum of ISO measurements were taken over the study site within 
100 m of the flux station during each mission. Two to five minutes were required to obtain these 
measurements. Data were collected under clear sky conditions one or two times a day. The sampled 
areas within the site were usually concentrated upwind of the micrometeorological flux sensors. The 
vegetation across the site was fairly uniform, with little variation in topography. 

The raw digital count values from the MMR were converted to sensor radiances using the 
calibration data taken prior to the experiment. The sensor radiances were then corrected for 
atmospheric effects using an atmospheric correction algorithm developed for FIFE by Fraser et al. 
(1989) based on the Dave (1978) algorithm that includes the effects of multiple scattering in the 
atmosphere and between the atmosphere and the surface. Aerosol optical depths were determined 
from a Sun-tracking photometer, operating in discrete bands in the 380 - 1020 nm region (Wriglcy 
et al., 1990 and Halthore et al. 1990) during the helicopter data acquisition. Gaseous absorption terms 
(CO* Oj, H 2 0, etc.) were calculated using Lowtran-7 (Kniezys, 1988) from radiosonde atmospheric 
water-vapor profiles taken simultaneously with the helicopter data acquisition. The SR was then 
calculated for the site using the average site reflectance in MMR3 and MMR4. 

2. 1.2.5 Estimating Canopy Conductance 

Canopy stomatal conductance (gj can be estimated from the surface conductance, g e (PM) 
[canopy conductance and surface conductance are roughly equivalent as the soil evaporation tends to 
zero - see e.g., Monteith, 1965; Thom, 1975; Baldocchi et al., 1991, Kim and Verma, 1991]. Values 
of gc(PM) = l/r c (PM) were calculated from the measured fluxes employing the Penman-Monteith 
equation (Monteith, 1965; Jarvis et al., 1976): 
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r c {PM) = [(s/y) p - 1] r a * 0+1) (pCJy) (D/A) , 


( 7 ) 


where s is the slope of the saturated vapor pressure-temperature curve, 7 is the psychrometric 
constant, p is the Bowen ratio (=HAE), H is the sensible heat flux, AE is the latent heat flux, p is the 
density of air, C p is the specific heat of air, D is the vapor pressure deficit, A (=R n +G) is the 
available energy, G is the soil heat flux and r a is the aerodynamic resistance to water vapor/heat 
transfer. The resistance, r, can be given by: 


r * = 


+ r L 


( 8 ) 


where is the aerodynamic resistance to momentum transfer (see Appendix B for details). The term 
r M can be expressed as: 

r am - ml , w 

where U is the mean horizontal wind speed at a reference height z and u* is the friction velocity. The 
excess resistance term, r„ (see e.g., Thom, 1972; Wesely and Hicks, 1977), is approximated by: 

« Df” . d») 

where k is von Karman’s constant, k is the thermal diffusivity and D v is the molecular diffusivity of 
water vapor. Friction velocity data from the three-dimensional sonic anemometer signals were used 
to compute r,. Values of gc(PM) were calculated using r, and measured fluxes of R n , G, H and AE in 
Eq. (7). 
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2. 1.2.6 Calculation of Canopy Photosynthesis (PJ1 

The eddy correlation measurement of atmospheric C0 2 flux (FJ provides the sum of canopy 
and soil C0 2 fluxes. An estimation of soil C0 2 flux is necessary for computation of canopy 
photosynthesis (PJ. Soil C0 2 flux was not directly measured in this study. However, our nocturnal 
F c data provided information on the respiratory release of CO z from plant and soil. Plant respiration 
(Rj) was estimated by the following relationship (Polley et al„ 1991): 

R d = R d (0 25° C) [exp iE t ( T t - 25)/(298 R r t )}] L T (11) 


where E, is the activation energy =45,000 J mol ', T, and T k are the leaf temperature in # C and °K, 
respectively, R is the gas constant, and Lp is the leaf area index. In our computations, the air 
temperature near the top of the canopy was used to approximate the leaf temperature. Following 
Polley et al. (1991), a value of 1.55 n mol m' 2 (leaf area) s' 1 was used for R d (@ 25°C). Plant 
respiration, estimated by Eq. (11), was subtracted from the nocturnal F c to obtain (nocturnal) soil C0 2 
flux. The soil C0 2 flux values were then adjusted to the daytime soil temperature (0-50 mm) using 
a Q l0 factor of 2 (Kucera and Kirkham, 1971) to calculate the daytime soil C0 2 flux. The daytime 
soil C0 2 flux was used in conjunction with the F c to estimate canopy photosynthesis (P c ). To minimize 
the effect of changing soil moisture and photosynthate on soil C0 2 flux, the nocturnal F c data used 
here were generally limited to those from the preceding or the following nights. 

2.1.3. Results and Discussion 

2.1.3.1 Soil Moisture and Plant Growth 

A seasonal picture of soil moisture and plant growth at the study site is presented in Fig. 1. 
Except for a three-week period in late July-early August, precipitation was generally ample during 
the study period. During the dry period, the extractable soil water, 6, (0-1.4 m), reached as low as 
25% and moisture stress conditions prevailed. Values of green leaf area index (LAI) and stages in the 
annual life cycle of the prairie vegetation are also included in Fig. 1. 
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2. 1.3.2 Seasonal Variations in Canopy Surface Conductance. Canopy Photosynthesis and 
Simple Ratio Vegetation Index 

In Fig. 2, an overall seasonal pattern 2 of canopy conductance, computed from our measured 
surface fluxes using Eqs. (7-10), is presented. During the most vigorous stage of plant growth (June 
and July), the midday g c averaged about 15.0 mm s' 1 . The dry period during late July-early August 
produced much lower value of g c with an average of 3 mm s' 1 . With ample rainfall in mid-August, 
the vegetation began to recover from moisture stress, resulting in an increase in g*. In mid-late 
August, the midday g c averaged about 9.0 mm s' 1 . Later in the season g, decreased as the prairie 
vegetation senesced. In October the midday g c averaged about 1.3 mm s' 1 . 

The seasonal variation 3 of canopy photosynthesis (P c ) is given in Fig. 3. The midday P e 
reached a seasonal peak of about 1.5 mg m' 2 (ground area) s' 1 during June-early July, the most 
vigorous stage of plant growth. There was a drastic reduction in P c during the dry period. The value 
of P c reached a minimum of near zero around 30 July, a day with severe moisture stress. The soil 
moisture conditions improved during mid-August, and the midday P e increased to about 0.75 mg m' 2 
s' 1 . The photosynthetic rates diminished rapidly as the vegetation senesced. The midday value of P c 
in early to mid-October was about 0.1 mg m' 2 s' 1 . 

The days on which SR was measured are indicated on Figs. 2 and 3. The seasonal progression 
of SR for the study site is shown in Fig. 4. Values of SR were high, ranging from 8.5 to 9.4, in June- 
early July when the plants were growing vigorously. During the early senescence stage which 
followed the dry period in late July-early August (see Fig. 1), the SR values reduced to 5. 1-6.2. The 
lowest values (~2) of SR were observed during mid-October toward the end of the senescence period. 

2.1.3.3 Simple Ratio Vegetation Index. Related to Unstressed Canopy Conductance and 
Unstressed Canopy Photosynthesis 

The procedure described in Appendix A was employed to define the form and coefficients 
for the stress factors, f(x), used in Eqs. 3 and 4. Measured values of D and 0, were used to calculate 


2 Details of these results have been reported in Kim and Verma, 1990a. 

3 Details of these results have been reported in Verma et al„ 1991. 
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the stress factors f(D) and f(0,). The values of g,. and P c were then divided by the stress factors to 
estimate g’ and P* (see Eqs. 3 and 4). Plots of g* and P’ against PAR are shown in Figs. S and 6. 
Estimates of dg^dPAR and dP^dPAR were obtained by applying simple linear fits to these data (only 
afternoon data were used in estimating BgJdPAR to minimize confounding effects of rapid time rate 
of change in soil evaporation). The values of dg^dPAR and dP^/dPAR are shown in Table 1. The 
correlation coefficient (r 2 ) was generally quite high, except on 4 June, 11 August and 8 October. 
Because of rain on 2 and 3 June, our analysis of data on 4 June is subject to problems arising from 
a rapid rate of change in soil evaporation. 11 August was a day with the driest soil moisture 
conditions of all the days during June-August when g* and P* were computed. Our empirical 
approach to calculate the stress factors under such severe moisture stress conditions is not adequate. 
The low r 2 on 8 October is likely because the canopy conductance and photosynthetic rates diminished 
rapidly as vegetation senesced in October. 

On five days (16 and 17 August and 7, 10 and 13 October), the canopy photosynthesis and 
conductance measurements were not available. The dg^dPAR and dP^/dPAR values on 16 and 17 
August were interpolated from the values on 15 and 20 August. On 7, 10 and 13 October, the values 
of dg'JdPAR and 3P^/dPAR were assumed to be the same as those on 8 October. 

In Figs. 7 and 8, SR is plotted against dg^/dPAR and dP^dPAR, respectively. As indicated 
above, the data included in these two figures are from a range of soil moisture and meteorological 
conditions and cover different stages in the annual life cycle of the prairie vegetation (Table 1). The 
data from 4 June-11 July are from a period when the soil water was not limiting and the plants were 
growing vigorously. The data point from 11 August is during the dry period, when severe moisture 
stress conditions prevailed. With about 67 mm of rainfall on 13 August, the soil water conditions 
improved significantly. The data collected during 15-20 August represent the improved soil water 
conditions. The data during 7-13 October are from the senescence period. 

Whereas the SR and canopy biophysical variables (dg^/dPAR, dP^/dPAR) track each other 
fairly well, the same cannot be said of the reported green LAI values and these three quantities. One 
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reason 4 may be that the canopy was much less lush in August than in June-July, but the bulk of the 
foliage was still classified as "green* rather than brown in August This would explain the lower SR 
values in August (in spite of larger values reported for green LAI) and also the lower values of 
dg^dPAR and dP'jdPAR. 

A linear regression of all of the SR-dP PAR data yielded a r 2 of 0.41 (Fig. 8). Similarly, 
a linear regression of all the data in the SR-dg ^/flPAR relationship (Fig. 7) resulted in r 2 =0.56. Our 
representation of the effect of moisture stress on g*and P’ is inadequate on 11 August when extreme 
dry soil moisture conditions prevailed. Also, as mentioned above, the value of dg^/dPAR on 4 June 
was probably affected by problems associated with a rapid change in soil evaporation. Excluding the 
data points on 4 June and 11 August a linear regression of the SR-dg ’Jd PAR data (solid line shown 
in Fig. 7) resulted in r 2 =0.69. Similarly, excluding the data point on 11 August, a linear regression 
of the SR-3P ’Jd PAR data (solid line shown in Fig. 8) yielded r 2 =0.57. 

Even though there is some scatter in the data. Figs. 7 and 8 seem to support the hypothesis 
suggested by Sellers (1987) that there should be near-linear relationships between the simple ratio 
vegetation index (SR) and the response of the unstressed canopy conductance and the photosynthetic 
capacity to the incident PAR (Eq. 5). The data scatter in Figs. 7 and 8 could have resulted from 
several simplifying assumptions used in the analysis presented above. First, there was no allowance 
for soil evaporation in the calculation of g c (PM); it was assumed that all of the latent heat flux was 
due to canopy transpiration. Second, the soil respiration was not measured directly, but was estimated 
using empirical relationships. Third, the expressions used to calculate the environmental stress factors 
are empirical and need further evaluation. 

The scatter in Figs. 7 and 8 is consistent with expectations given the experimental 
methodologies and simple analytical tools we have applied. Further work is needed to address the 
problems outlined above to bring us closer to the goal of using satellite data to provide information 
on energy and mass exchanges in various terrestrial ecosystems. 

*To estimate LAI, leaves of eight 0.18m 2 area were harvested. The harvested foliage was divided 
into live and dead components, and then fed through an area meter. However, this procedure took 
no account of changes in the greenness/lushness of live vegetation and in August the leaves were less 
green than in June/July. 
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2.1.4. Summary and Concluding Remarks 

The surface fluxes of carbon dioxide and water vapor were measured employing the 
micrometeorological eddy correlation technique during FIFE in 1987. Concurrent measurements were 
also made of reflectance using a helicopter-mounted Barnes MMR, from which the simple ratio 
vegetation index (SR) was calculated. Both of these data sets provide information on the grassland 
canopy at a comparable scale, (10-100 m) 2 . 

Using our data, both on the leaf level (e.g., stomatal conductance of the three dominant 
species) and the canopy level (e.g., water vapor flux), empirical relationships for the environmental 
stress factors, f(D) and f(0,), were developed to account for the effect of ambient vapor pressure 
deficit (D) and extractable soil water (8,). These factors were then used in conjunction with the 
measured surface fluxes to compute unstressed canopy conductance (g*) and unstressed canopy 
photosynthesis (P’). 

The derivatives (dg^/3PAR and dP^/dPAR) of unstressed canopy conductance and unstressed 
canopy photosynthesis are dependent on the amount and type of green vegetation, and therefore, 
should be amenable to remote sensing. The analysis of Sellers (1987) indicated that dg^/dPAR and 
dP^/dPAR should be near-linearly related to SR and the first order, preliminaty analysis of data 
presented in this paper seems to support this hypothesis. However, before such relationships can be 
employed in satellite remote sensing applications, more work is needed. Future studies should 
include: (a) adequate allowance of soil evaporation in calculating canopy conductance, (b) direct 
measurement of soil respiration, (c) further development and evaluation of environmental stress 
functions and (d) more frequent measurements of canopy reflectance under similar illumination 
conditions. 

2.2 Further Data Processin g and Analyses 

All of the eddy correlation data collected at site 916 during 1989 have been processed and 
fluxes of energy and mass have been calculated. Supporting meteorological, soil water and plant 
growth data have also been computed. 
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Soil C0 2 flux has been estimated using the relationships developed at our site by Norman el 
al. (1991). Employing these estimates in conjunction with our eddy correlation measurements during 
FIFE-89, canopy photosynthesis of the prairie vegetation is being estimated. A biochemical model 
of leaf photosynthesis was applied to the tallgrass prairie. The modeled photosynthetic rates of 
individual leaves were scaled up to the canopy level using a radiative transfer model. The FIFE-89 
data are now being used to further test the canopy photosynthesis model. This model as well as a 
’multiplicative" stomatal conductance model are employed to evaluate canopy conductance. Further 
tests of these models using the FIFE-89 data are underway. 

3. JOINT STUDIES WITH OTHER FIFE INVESTIGATORS 

The following includes a list of major joint studies presently underway: 

• Relating surface fluxes to reflectance values measured by the helicopter-mounted 
multiband radiometer (MMR) — with Dr. G Walthall, University of Maryland, College 
Park, MD, Drs. P. Sellers and F. Hall, NASA-Goddard Space Flight Center, Greenbelt, 
MD. 

• Comparison of canopy surface conductance at two contrasting sites within the FIFE 
study area — with Dr. J. B. Stewart, Institute of Hydrology, Wallingford, Oxon, United 
Kingdom. 

• Intercomparison of rapid response C0 2 sensors — with Dr. M. Wesely, Argonne National 
Laboratory, Argonne, IL and Dr. J. Moncrieff, University of Edinburgh, Edinburgh, 
U.K 

4. SCIENTIFIC PUBLICATIONS AND PRESENTATIONS STEMMING FROM GRANT 
NAG5-890 

The following scientific publications and presentations stem from our work on Grant No. NAG5-890. 

Journal Articles 

Hall, F. G., P. J. Sellers, L McPherson, R. D. Kelly, S. B. Verma, R. Markham, B. L. Blad, J. Wang 
and D. E. Strebel. 1989. FIFE: Analysis and Results — A Review. Advances in Space 

Research. Cospar Publication 9:275-293. 
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Kim, J. and S. B. Verma. 1990. Components of surface energy balance in a temperate grassland 
ecosystem. Boundary- Layer Meteorology 51:401-417. 

Kim, J. and S. B. Verma. 1990. Carbon dioxide exchange in a temperate grassland ecosystem. 
Boundary-Layer Meteorology 52:135-149. 

Verma, S. B. 1990. Micrometeorological methods for measuring surface fluxes of mass and energy. 
In: Instrumentation for Studying Vegetation Canopies for Remote Sensing in Optical and 
Thermal Infrared Regions, N. S. Goel and J. M. Norman (co-editors). Special Issue of Remote 
Sensing Reviews 5(1):99-1 15. 

Kim, J. and S. B. Verma. 1991. Modeling canopy stomatal conductance in a temperate grassland 
ecosystem. Agric. and Forest Meteorol. 55:149-166. 

Kim, J. and S. B. Verma. 1991. Modeling canopy photosynthesis: Scaling up from a leaf to canopy 
in a temperate grassland ecosystem. Agric. and Forest Meteorol. (accepted for publication). 

Verma, S. B., J. Kim and R. J. Clement. 1991. Momentum, water vapor and carbon dioxide 

exchange at a centrally located site in FIFE. J. Ceophys. Res. (submitted). 

Stewart, J. B. and S. B. Verma. 1991. Comparison of surface fluxes and conductances at two 
contrasting sites within FIFE area. J. Geophys. Res. (submitted). 

Moncrieff, J. B., S. B. Verma and D. R. Cook. 1991. Intercomparison of eddy correlation carbon 
dioxide sensors during FIFE-89. J. Geophys. Res. (submitted). 

Norman, J. M., R. L. Garcia and S. B. Verma. 1991. Soil C0 2 fluxes and the carbon budget of a 
grassland. J. Geophys. Res. (submitted). 

Smith, E. A., A. Y. Hsu, W. L. Crosson, R. T. Field, L. J. Fritschen, R. J. Gurney, E. T. Kanemasu, 
W. P. Kustas, D. Nie, W. J. Shuttleworth, J. B. Stewart, S. B. Verma, H. L. Weaver and M. 
L. Wesely. 1991. Area averaged surface fluxes and their time-space variability over the 
FIFE experimental domain. J. Geophys. Res. (submitted). 

Conference Presentations 

Verma, S. B. 1988. Carbon dioxide and water vapor fluxes over a tallgrass prairie during FIFE. 
Invited Paper. American Geophysical Union Spring Meeting, Baltimore, MD, May 18-20, 
1988. 

Sellers, P. J., F. G. Hall, D. E. Strebel, R. D. Kelly, S. B. Verma, B. L. Markham, B. L. Blad, D. S. 
Schimel, J. R. Wang, W. Brutsaert, R. Desjardins, L. Fritschen, R. Frouin, R. Grossman, E. 
Kanemasu, L McPherson, J. L. Dorman and G Walthall 1988. First ISLSCP Field 
Experiment: Experiment Execution and Preliminary Analyses. European Center for Medium 
Range Weather Forecasts (ECMWF) Symposium on land surface parameterizations, 
Proceedings, October 1988. Shinfield, U.K. 

Kim, J. and S. B. Verma. 1989. Energy and C0 2 Fluxes in a Grassland Ecosystem. Presented at the 
19th Conference on Agricultural and Forest Meteorology of the American Meteorological 
Society, Charleston, SC, March 7-10, 1989. 

Norman, J. M., W. Polley, E. A. Waiter-Shea, S. B. Verma, T. Arkebauer and D. Greegor. 1989. 
Comparison of measured and modeled fluxes of water and carbon dioxide at Konza Prairie. 
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Presented at the 19th Conference on Agricultural and Forest Meteorology of the American 
Meteorological Society, Charleston, SC, March 7-10, 1989. 

Stewart, J. B. and S. B. Vertna. 1989. Comparison of aerodynamic and surface conductance at two 
contrasting sites within FIFE area. Presented at the Global Energy and Water Fluxes 
Symposium, Fifth Scientific Assembly of the International Association of Meteorology and 
Atmospheric Physics (IAMAP), Reading, U.K., July 31 -August 12, 1989. 

Kanemasu, E T., S. B. Verma, L J. Fritschen and R. J. Gurney. 1990. Design and operation: flux 
measurements in FIFE Presented at the American Meteorological Society Symposium on the 
First ISLSCP Field Experiment (FIFE), 70th AMS Annual Meeting, February 7-9, 1990, 
Anaheim, CA. 

Norman, J. M., S. B. Verma and R. Garcia. 1990. Soil surface C0 2 fluxes in the Konza Prairie. 
Presented at the American Meteorological Society Symposium on the First ISLSCP Field 
Experiment (FIFE), 70th AMS Annual Meeting, February 7-9, 1990, Anaheim, CA. 

Stewart, J. B. and S. B. Verma. 1990. Comparison of surface conductance at two contrasting sites 
within FIFE study area. Presented at the American Meteorological Society Symposium on 
the First ISLSCP Field Experiment (FIFE), 70th AMS Annual Meeting, February 7-9, 1990, 
Anaheim, CA. 

Verma, S. B., J. Kim and R. J. Clement. 1990. Fluxes of momentum, sensible heat, water vapor and 
C0 2 over a tallgrass prairie. Presented at the American Meteorological Society Symposium 
on the First ISLSCP Field Experiment (FIFE), 70th AMS Annual Meeting, February 7-9, 
1990, Anaheim, CA. 

Stewart, J. B. and S. B. Verma. 1991. Can the simple ratio be used to estimate evaporation from the 
FIFE area? To be presented at the Special Session on Hydrometeorology of the American 
Meteorological Society, September 10-13, 1991, Salt Lake City, UT. 

Kim, J. and S. B. Verma. 1991. Modeling canopy photosynthesis in a temperate grassland ecosystem. 
To be presented at the 20th Conference on Agricultural and Forest Meteorology of the 
American Meteorological Society, September 10-13, 1991, Salt Lake City, UT. 


5. RESEARCH PLANS 

Estimation of canopy photosynthesis during the measurement period in FIFE-89 will be 
completed. Further testing of canopy photosynthesis model, presently underway, will be completed 
in the next few months. Results from the 1989 stomatal conductance measurements will be 
incorporated into a multiplicative canopy conductance model. This model will also be further tested 
with FIFE-89 data. 

Our micrometeorological measurements during FIFE-87 and FIFE-89 will be used to 
calculate daily net ecosystem C0 2 exchange. These results will be used to examine the contribution 
of this temperate grassland as a source and sink for atmospheric C0 2 under a wide range of 
meteorological and soil moisture conditions. 
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7. 


APPENDICES 


APPENDIX A 


The procedure employed to derive the stress factors f(T), f(D) and f(0J (see Eqs. 3 and 4) is 
described here. 

Leaf stomata! conductance (g,) can be described in terms of a simple multiplicative model (see 
e.g., Jarvis, 1976; Kim and Verma, 1991): 

g, - g,<PAK> ■ m m wi. < A -» 


where g,*(PAR) is the unstressed value of leaf stomatal conductance. Similarly, leaf photosynthesis 
can be estimated as: 


p t = p t \PAR) • \AD m y(0;] , 


(A-2) 


where P,*(PAR) is the unstressed value of leaf photosynthesis. The stress factors f(T), f(D) and f(tf,) 
account for the effect of temperature (T in C), ambient vapor pressure deficit (D in kPa) and 
extractable soil water (8, in %). Values of were computed as the ratio of actual to total soil moisture 
held with a water potential between -1/30 and -1.5 MPa. 

Equations (A-l) and (A-2) can be integrated ("scaled up") to provide estimates of canopy 
stomatal conductance (gj and canopy photosynthesis (P c ); 

g c = F(E) f h g ;(PAK)dL = g c '(PAK) • F(E) (A-3) 

» O 


and 


21 



P c = F(E) j L J P,'(PAK)dL = P C '(PAK) • F(E) , 


(A-4) 


where F(E) = f(T)«f(D)*f(0J and Lj. = total leaf area index. We assume that the stress factors are 
independent of leaf area index. 

There appears to be a wide range of temperature optima for the grass species in the 
experimental area (Knapp, 1985) and, therefore, the effect of changes in T was assumed to be 
negligible (i.e., f(T) » 1). The response of the unstressed leaf stomatal conductance to PAR was 
estimated using a hyperbolic relationship (e.g., Monteith, 1965): 

8,* (PAR) = (a^AR i)/(a 2 + PAR 1) , (A-5) 


where PARI is the normal (relative to the leaf) flux density of incident PAR, and a, and a 2 are 
empirically derived constants 5 . The leaf stomatal conductance data of the three dominant species (big 
bluestem, indiangrass and switchgrass) measured at the site were used to evaluate parameters a, and 
a 2 for each of the three species (see Kim and Verma, 1991 for details). Averaged values of these 
parameters (a 2 = 19.7 mm s 1 and a 2 = 441 /rEi m' 2 s ') were used here to yield a canopy level 
relationship for the mixed vegetation. 

Using f(T) = 1 and Eq. (A-5), Eq. (A-3) can be rewritten as: 


8 C 



r a { PARl ' 


dL 

02 + PAR\ t 


AD) AV • 


(A-6) 


Substituting the Ross-Goudriaan function (see e.g.. Sellers, 1985, Table 3, p. 1351), Eq. (A-6) can be 
simplified to: 


5 The parameter, a„ represents the asymptotic value of g ( * when PAR-«°, and can be considered 
as the maximum stomatal conductance under optimal conditions; and a 2 defines the curvature of the 
response to PAR. 
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Sc * 


tot- 


Hflj + G(n) PAR 0 


V<h + PAR o e 


-*,b 


7(D) /(e; , 


(A-7) 


where k„ = extinction coefficient for PAR =(1 - w,) w G(ji)lp, w, - 0.2, G (ji) «0.5 (for spherical leaf 
distribution) and PAR„ = PAR flux above the canopy. 

A curvilinear reduction in canopy conductance with increasing vapor pressure deficit (D) was 
assumed (e.g., Lohammar et al„ 1980): 

/(D) = 1/(1 + a } D) , (A- 8 ) 


where a 3 is a constant. The response of canopy conductance to 8, (extractable soil water) was 
estimated using a negative exponential relationship (e.g., Stewart, 1988; Kim and Verma, 1991): 

w - • - . <*-») 


where a 4 is a constant. 

As discussed in Sec. 2.5, g,. was computed from our half-hourly measurements of fluxes of net 
radiation, soil heat, sensible heat and latent heat. Values of n and k, were calculated every half hour. 
The stress functions f(D) and ((8,) were then evaluated by fitting (nonlinear least squares) these data 
(&, *i. k„ PAR e , L,.) to Eq. A-7: 

03 = 0.2031 ± 0.0536 

0 4 = 0.0144 ± 0.0017 (A- 10) 

(r 2 = 0.788) 
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(These values were derived from 4SO half-hour & data points collected during 20 days distributed 
throughout the season.) 

Thus, the following stress factors were used in this study: 

f(D) = 1/(1 +0.2031D) < A - n ) 


and 


= 1 - «p(-0.0144e,) 


(A-12) 


where D is the air vapor pressure deficit measured at 2.25 m above ground and $, is the daily value 
of extractable soil water over the primary root zone (0 to 1.4 m depth). 

These derived stress factors were applied to the estimates of g e and computed from flux 
observations, to calculate values of g’ and P* (see Sec. 3.3). Note that while values of L,. were used 
to derive the coefficients given in Eqs. (A-ll) and (A-12), the coefficients and the forms of stress 
factors are independent of the time-varying leaf area index. 

APPENDIX B 


Heat (or water vapor) transfer in the immediate vicinity of a vegetated surface is primarily 
controlled by molecular diffusion. Momentum, however, is transferred to a large extent by pressure 
forces, which have no analog in heat (or water vapor) transfer. Transfer of momentum in 
aerodynamically rough flows (e.g., over a vegetated surface) is primarily independent of molecular 
viscosity. Therefore, the transfer of heat (or water vapor) to or from a vegetated surface encounters 
greater aerodynamic resistance than does the transfer of momentum (Owen and Thomson, 1963; 
Chamberlain, 1966, 1968; Thom, 1972; Verma, 1989). Accordingly, an 'excess resistance’ term, r b is 
defined: 
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where r^ and r, are the aerodynamic resistance to transfer of momentum and sensible heat/water 


vapor. 

The terms r,^ and r a can be obtained from: 



and 


H = p 



(B-2) 


(B-3) 


where r (=p u* 2 ) is the momentum flux, u* is the friction velocity, 0 and T are the mean horizontal 
wind speed and air temperature at a reference height The term T a represents an effective surface 
temperature. Equation (B-2) can be rewritten as: 



(B-4) 


Owen and Thomson (1963) proposed a dimensionless parameter B' 1 (sublayer Stanton number) 
to express the difference between r, and r^,: 

B l = r* - r* , (B-5) 

" a 'am » 


where r/ and r ao + are nondimensionless resistances given by: 


and 



(B-7) 


Several investigations (e.g., Garratt and Hicks, 1973; Brutsaert, 1984) suggest that for surfaces with 
densely spaced permeable roughness elements, as a first approximation k B' 1 can be assumed to be of 
the order of 2 for scalars whose Schmidt (S c ) or Prandtl (P f ) number is of the order of 0.6 to 0.8 (k is 
von Karman’s constant). Therefore, the excess resistance term, r* can be expressed as: 

r h « (2/fcO [SJPy* = (2/te.) (k JDf 3 . (B-8) 


The factor (S,/P r ) w is incorporated to account for the fact that the basic information was derived from 
heat transfer observations primarily (Wesely and Hicks, 1977). The term k represents the thermal 
diffusivity and D v represents the molecular diffusivity of water vapor. 


26 



APPENDIX C 


Table 1. Values of dg^/dPAR and AP^/dPAR. 


Date (1987) 

g\ vs PAR 

p; ps PAR 1 

(dg^dPAR)xlOOO 
(mm m'V*) 

r* 

(dP^dPAR)xlQ00 
(mg m-VVEI nr*-*) 

r* 

June 4 

14.0 

0.590 

1.01 

0.845 

June 5 

9.4 

0.976 

1.23 

0.951 

June 6 

8.6 

0.964 

1.19 

0.927 

July 11 

10.8 

0.887 

1.24 

0.966 

August 11 

0.2 

0.002 

0.14 

0.029 

August IS 

8.2 

0.961 

1.16 

0.961 

August 16* 

8.6 

— 

1.18 

— 

August 17* 

9.1 

— 

1.19 

— 

August 20 

10.3 

0.955 

1.24 

0.923 

October 7* 

2.5 

— 

0.44 

— 

October 8 

2.5 

0.438 

0.44 

0.683 

October 11* 

2.5 

... 

0.44 

— 

October 13* 

2.5 

... 

0.44 

... 


'Estimated data (see text for details). 
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Table 2. Information of annual life cycle stage, green leaf area index, extractable soil water ( 9 
0-1.4 m depth), midday (averaged during 1230-1430 hrs, central daylight savings lime) 
mean air temperature (T) and vapor pressure deficit (D) and photosynthetically active 
radiation (PAR) on selected days when reflectance and flux measurements were 
available. T and D were measured at 2.25 m above ground. 


Date 

Day 

Number 

Annual Life 
Cycle Stage 

E9 

Green 

LAI 

H 

D 

(KPa) 

PAR 
OiEl in'* 
•- 1 ) 

June 04 

155 

Early Growth 

79.3 

1.86 

25.4 

1.88 

2052 

June 05 

156 

Early Growth 

77.7 

1.90 

27.1 

1.91 

2048 I 

June 06 

157 

Early Growth 

76.2 

1.94 

27.8 

2.11 

2050 | 

July 11 

192 

Peak Growth 

68.4 

179 

31.0 

1.85 

1971 

Aug. 11 

223 

Early 

Senescence 

24.8 

2.44 

32.6 

2.87 

1844 

Aug. 15 

227 

Early 

Senescence 

54.4 

2.39 

33.7 

2.35 

1888 | 

Aug. 16 

228 

Early 

Senescence 

54.2 

2.38 

30.8 

1.57 

1887 1 

Aug. 17 

229 

Early 

Senescence 

54.0 

2.37 

31.4 

2.71 

1885 I 

Aug. 20 

232 

Early 

Senescence 

53.4 

2.28 

32.8 

2.35 

1860 | 

Oct. 07 

280 

Senescence 

30.7 

0.30 

15.9 

1.71 

1351 | 

Oct. 11 

284 

Senescence 

29.8 

0.20 

11.6 

0.97 

1450 | 

Oct. 13 

286 

Senescence 

32.5 

<0.20 

22.7 

2.13 

1165 | 
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APPENDIX D 


LIST OF FIGURES 

Fig. 1. Seasonal precipitation, extractable soil water (0-1.4 in depth) and green leaf area index (after 
Kim and Verma, 1991). 

Fig. 2. Seasonal variation in midday (1230-1430 h, central daylight savings time) values of canopy 
conductance, (after Kim and Verma, 1990a). The days on which SR was measured are 
indicated with arrows. 

Fig. 3. Seasonal variation in midday (1230-1430 h, central daylight savings time) values of canopy 
photosynthesis (after Verma el al., 1991). The days on which SR was measured are indicated 
with arrows. 

Fig. 4. Seasonal variation in the simple ratio vegetation index (SR). 

Fig. S. Unstressed canopy conductance (g*) as a function of the photosynthetically active radiation 
(PAR). Only afternoon data were used to minimize confounding effects. Estimates of 
dgJdPAR were obtained by applying simple linear fits to these data. 

Fig. 6. Unstressed canopy photosynthesis (P*) as a function of the photosynthetically active 

radiation (PAR). Estimates of 3P‘ /a PAR were obtained by applying simple linear fits to 
the data. 

Fig. 7. Simple ratio vegetative index (SR) as a function of the derivative (dg^/dPAR) of the 

unstressed canopy conductance with respect to the photosynthetically active radiation (PAR). 
The correlation coefficient (r 2 ) for the linear regression is 0.56 for all data and 0.69 for the 
data set excluding 4 June and 11 August. The solid line represents the linear regression of 
the data set excluding 4 June and 11 August. 

Fig. 8. Simple ratio vegetative index (SR) as a function of the derivative (aP^dPAR) of the 

unstressed canopy photosynthesis with respect to the photosynthetically active radiation 

(PAR). The correlation coefficient (r 2 ) for the linear regression is 0.41 for all data and 0.57 
for the data set excluding 11 August. 
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